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Abstract 

The Busch-formula relates the energy-spectrum of two point-hke particles interacting in a 3-D isotropic 
Harmonic Oscillator trap to the free scattering phase-shifts of the particles. This formula is used to find an 
expression for the shift in the spectrum from the unperturbed (non-interacting) spectrum rather than the 
spectrum itself. This shift is shown to be approximately A = —5{k)/-n x dE, where dE is the spacing between 
unperturbed energy levels. The resulting difference from the Busch-formula is typically < \% except for the 
lowest energy-state and small scattering length when it is 3%. It goes to zero when the scattering length 
— ±cx). 

The energy shift A is familiar from a related problem, that of two particles in a spherical infinite square- 
well trap of radius R in the limit R ^ oo. The approximation is however as large as 30% for finite values of 
7?, a situation quite different from the Harmonic Oscillator case. 

The square-well results for i? — >■ oo led to the use of in-medium (effective) interactions in nuclear matter 
. calculations that were oc A and known as the ptiase stiift approximation . Our results indicate that the 

04 ' validity of this approximation depends on the trap itself, a problem already discussed by DeWitt more than 

50 years ago for a cubical vs spherical trap. 
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1 Introduction 



In some of his early attempts of developping a many-body theory of nuclei, Brueckner assumed the in-medium 
(effective) two-body interaction to be cx tanS{k), S{k) being the scattering phase-shift as a function of relative 
momentum k. The initial problem studied was that of an 'infinite' system of nuclear matter for the purpose of 
calculating binding energy and saturation properties. The tan(5(fc) approximation came from the assumption 
O ' that the Reactance matrix, a part of scattering theory, could be used as an in-medium interaction in the many- 

2 . body problem. (The diagonal part of this matrix is cx ta,ii6(k)). This implies that a principal value Green's 

function with an integration over a semi-continuous spectrum would be justified in this case with the box of 
nuclear matter assumed to be very and even 'infinitely' large. This assumption was substantiated in a paper by 
Reifman and DeWitt[T]. It was however soon realised to be incorrect. Several authors (one of them DeWitt) 
showed that the correct limit for two particles in a big box would yield an in-medium interaction oc S{k) rather 
than tan 5(A:). [3 H S] 1 The Reactance matrix is part of scattering theory. The nuclear matter problem 
. assumes a box, although large but still finite with boundary conditions different from that of scattering theory. 

' A particle in this box has a discrete rather than a continuous spectrum and that makes a difference.. 

The proofs in the referred papers differ but the essential point is that the integration over the continuous 
spectrum vs the summation over the discrete spectrum differs even though the level-density in a big box goes to 
zero as the box-size increases. The principal value integration relates to the scattering problem with boundary 
conditions different from that for a box where the wave-functions are zero at the edge of the box even in the 
limit of 'infinitely' large. 

The exact statement of the results in the referenced papers is that the energy-shift A due to the interaction 
^ j , of two-particles confined in a spherical box in the limit when the size of the box approaches infinity is given by 

d ' A = S X dE where dE is the spacing between the levels of the unperturbed spectrum and H 

s^JM. (1) 

TT 

In nuclear matter studies this has been referred to as the 'phase-shift approximation' a term adopted here. It 
was for example used in early calculations on the neutron-gas [5l|6] and used as a first approximation for nuclear 
matter studies. [T 

The question of the author arose whether the result above might be more general. The problem of two 
particles trapped in a potential well is of interest for atomic as well as nuclear physics studies. The Busch- 
formulaiSj relates the energy- spectra of the two particles in a Harmonic Oscillator well to the scattering phase- 
shifts for point-like potentials. This formula is here rewritten in terms of the energy- stiifts rather than energy. 

^e-mail: kohlersCSu. arizona.edu 

^For small S e.g. at low density this might not make any difference but for problems of interest to-day with 5 7r/2 it obviously 
would. 

^It was pointed out by DeWitt[3], after a comment by Brueckner that the result may be different for a box other than spherical. 
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Results below show that the phase-shift approximation for these shifts is (surprisingly) in practically exact 
agreement with the Busch-formula. The largest difference is found for the lowest energy-state {n = 0) and 
small scattering length but it decreases rapidly for n > and larger scattering lengths. For comparison are 
also shown results for the spherical infinite square well as a function of the radius of the sphere. Although the 
approximation becomes exact in the limit R — !■ oo, it is in general much worse in this case than for the H.O. 



2 3-D Harmonic Oscillator Well 

With the energy in units of hw the total energy for two non-interacting particles with zero angular momentum 
in a 3-D Harmonic Oscillator well is 

Etot ~ E + Ecm = 2n + 3 
where Ecm = § is the center-of mass energy . It is convenient to choose 



as the unit of length. With rj = 2E the Busch-formula reads [3] 

tanJffc)-^-^^™^^ (2) 
taud[k)- 2r[(3-,7)/4] 

where k — y^/aosc- 

The level-spacing in the uncorrelated system is dE = 2 and following the notation above one finds the 
level-shift to be A = 25*. In the case of the spherical box above S is given by eq. ([Ij. It is the purpose of the 
present work to find S for the 3-D H.O. with the Busch spectrum. One finds 

77 = 4n + 3 + 45* 

. Using the reflection formula for F-functions 

F(l - x) = , J , , (3) 

^ ^ F(x) sin(7r2:) ^ ' 

one finds 

tan (5(fc) = -A{z) tan(S'7r) (4) 



where 



with z=n-|-S-|-l. 

The beta-function B{x,y) is defined by 



and for y fixed and x large one has 0. 
so that for z large one would have 



Mz) = V, ^1, (5) 
F(z+ 2) 



Bi.,y)^^J^ (6) 



T{x + y) 



Bix,y)^Tiy)x-y (7) 



r(z + h] 



F(z) 



but in cq. ([S]) one has the factor 
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''I owe my thanks to Dr Jerry Yang for showing me this relation on Wikipedia 
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Figure 1: The function A{z) defined in the text. 



Numerical tests show that a better approximation indeed is 



■ - M Z 

r(z) V 4 

One finds that the values of z that is needed in the Busch-formula is z > i and Fig. [1] shows the function 
A{z) for a useful range of z. 

One finds ^ — >■ 1 quite rapidly with z (and consequently with n) so that in this limit cq. (|4]) yields 

S^-'-^ (8) 

TT 

which is the phase-shift approximation given in eq. ([1]). 
A correction dS to S due to A{z) 7^ 1 is obtained from 

1 / , 1 / NN tan^(fc) ,„. 
'^^"-^(^-^ (-» l + tan^.(fc) 

This shows a dampening of the correction with increasing value of tan^ 6. In the limit when <5 = f the correction 
goes to zero. The energy-shift is then exactly given by = i, a well-known result.^ Fig. [5] shows the 
difference between the energy shifts calculated by the Busch-formula and the phase-shift-approximation for 
quantum levels n = 0, 1 and 2 as a function of scattering length. It is seen to decrease rapidly with n and also 
with ±a. The peaks are a result of the competition between the A[z) and tan (5 corrections respectively shown 
in eq. ([9]). Fig. [3] shows the same differences but in % of the Busch result for n = 1,2,3 again showing the 
rapid decrease with increasing values of n and with ±a. 



3 Spherical Square Well 

The problem of two particles in a spherical box was treated by several authors referred to in the Introduction. 
In these several works the emphasis was on the energy-shift in the limit of the radius of the sphere going to 
'infinity' with the result given by eq. ([I}, the phase-shift approximation. For comparison with the 3-D H.O. it 
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Figure 2: The difference between the energy-shifts calculated with the Busch- formula and that obtained in the 
phase-shift approximation for n = 0, 1, 2. The largest difference is as expected for n = (red on-line), while it 
is appreciably smaller for n — 1 (green on-line) and n = 2 (blue on-line). The smallest value of z in the function 
A{z) is i (see Fig. [1]), which occurs for n = and S — The largest difference shown is consequently for 
this value of n and for a < where the shift is negative. 




Figure 3: Similar to Fig. [2]but the difference is here shown as the percentage of the energy-shift. The maximum 
difference is « —3% (for n — 0, red on-line) but a rapid decrease with increasing n is seen and is less than i% 
for n — 1 (green on-line) and even smaller for n = 2 (blue on-line). 
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is of interest to also display the shift as a function of the ratio of scattering length to radius and as a function 
of quantum-numbers n of the spherical box. 

The simplest solution [21 [10] of the problem at hand is to explicitly consider the wave-functions of the two 
particles in this box. Considering only s-states the radial wavefunctions of free non-interacting particles are 



^-fr) oc isin(fc("V) (10) 
r 



With two-body interaction the wave-function outside the range of the two-body interaction is 

^'(r') cx - sin(fcr -|- 5{k)) (11) 

The boundary condition implies that the wave-functions vanish at the boundary of the sphere assumed to 
have a radius R. (One immediately sees the simplicity of the formalism relative to that of the 3-D H.O. or a 
cubical box.) Thus, for the non-interacting case 

= rnr (12) 

and for the interacting case 

knR + S{kn) = mr (13) 

The spacing between unperturbed levels is 

.(0) 2 



The energy shift 



^-kl-k(f^-{k,, + ki^^) 
divided by the unperturbed energyspacing, i.e. S, is then 

_ kn + kn / 5{kn ) \ 

~ + \ ^ J 

The energy-shifts and -spacings decrease with R and n ^ oo and then 

TT 

as before. 

The phase-shift approximation was tested for convergence as a function of box radius R and quantum number 
n. For this purpose the expression for S was rewritten as follows 

^_ 1 / S{k^)-2n7T) \ 6{k^) 

TT V 2n y TT ^ ' 

The radius R only appears implicitly through the relation 

kn = (nTT - 5{kn))/R 

to be solved selfconsistently. 

In the (unphysical) case that 5 would be independent of k one would have the situation where S is independent 
of i?, but only of n and of course of 5. 

Fig. [4| shows the difference between the 'exact' shift and the phase-shift approximation — |: of eq. ([T]), while 
Fig. [Sj shows the difference in % of the exact shift. 

Compared with the analogous result for the H.O. there are notable differences. While the largest deviation 
from the phase-shift approximation is 3% (and typically much smaller) in the H.O. case, it can for the sphere 
be even as large as 30%. But even though the devation from the phase-shift approximation in the H.O. is small 
it has its maximum for large values of aosc (or small scattering lengths) while going to zero for small values (or 
large scattering lengths), which is opposite to the situation encoutered for the spherical well. Both Figs. ID and 
[5| illustrate however the known result that S ^ — as i? — > oo as shown in previous works already referred 
to above, d [311]. 
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Figure 4: The difference between the exact and phase-shift approximation for two-particles in the spherical trap 
for three values of the quantum- number n — 1, (red on-line), n — 2 (green on-line) and n — 3 (blue on-line). 
Notice the difference with the similar plot in Fig. [2] for the H.O. trap. 




Figure 5: Similar to Fig. |4]but the difference is shown as percentage of the energy-shift. 
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4 Summary 



The phase-shift approximation —^d{k)dE for the energy-shift due to the interaction of two partcles in a trap, 
where d{k) is the free particle scattering phase-shift and dE the energy-level spacing of the unperturbed spec- 
trum, has been investigated. 

The traps considered were a 3-D Harmonic Oscillator and a spherical infinite square well. The square-well 
case was treated a long time ago as referenced above, with the main emphasis on the large size limit of the trap 
with the result expressed by eq. ([T]). Somewhat of a surprise was the finding of the present investigation that 
a not only somewhat similar situation exists for the H.O. trap but that the approximation is so much more 
accurate for this trap. 

The differences between the results for the two respective traps are seen by comparing Figs [2] and [3] for 
the H.O. trap with the Figs. |4] and [5] for the square-well trap. While in the square well case, the differences 
(errors) by using the phase-shift approximation are as large as 30%, the differences are much smaller and even 
practically zero in the H.O. case. With a — >■ ±oo the phase-shift approximation becomes exact in the H.O. case 
while in the spherical case the error increases. While Figs |4] and [5] show that the approximation becomes perfect 
when a/R goes to zero (or the radius of the square well goes to oo) that is when (unexpectedly) the largest 
correction, although still small, is found in the H.O. case. 

In all the work above like in the Busch-formula, it is assumed that the range of the potential is small 
compared to the size of the trap. Considerations of range corrections are beyond the scope of this investigation. 
It has been the subject of a recent publication. [TT| 

The question arises of course whether the phase-shift approximation has a more general validity e.g. for the 
shifts in a cubical box, as derived by Liischer.[12^ It is however already seen that the situation is rather different 
in the 3-D H.O. case compared to the case of a spherical well. The conceptual differences between scatterings 
and interactions in the spherical vs cubical box was already pointed out by DeWitt in a detailed discussion of 
this matter. [3]. 
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^ Quote: Only the spherical box (with spherical waves) is suitable for establishing a connection between single scattering processes 
and discrete spectrum theory.... 
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